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Abstract 

We compute, based on density-functional electronic-structure calculations, the Coulomb 
couplings in the hu highest occupied orbital of molecular Ceo • We obtain a multiplet-averaged 
Hubbard {7 ~ 3 eV, and four Hund-rule-like intra-molecular multiplet-splitting terms, each of 
the order of few hundreds of me Vs. According to these couplings, all Cg^ ions should possess a 
high-spin ground state if kept in their rigid, undistorted form. Even after molecular distortions 
are allowed, however, the Coulomb terms still appear to be somewhat stronger than the pre- 
viously calculated Jahn- Teller couplings, the latter favoring low-spin states. Thus for example 
in CgJ, unlike C^q , the balance between Hund rule and Jahn Teller yields, even if marginally, 
a high-spin ground state. That seems surprising in view of reports of superconductivity in 
field-doped C^q systems. 

1 Introduction 

Strong electron correlations in multi-band, orbitally degenerate systems represent an important 
current theoretical challenge. A lively experimental playground for that is provided by electron- 
doped fullerene systems, which exhibit a variety of behavior, including unconventional metals like 
cubic CsCeo (Brouet et al. 1999), superconductors of the A3C60 family (A= K,Rb,Cs) (Ramirez 
1994, Gunnarsson 1997), and insulators, presumably of the Mott- Jahn- Teller type (Fabrizio et al. 
1997, Capone et al. 2000), such as Na2C6o (Brouet et al. 2001), A4C60 (Benning et al. 1993), and 
the class of ammoniated compounds (NH3)K3_^Rbi,C6o. 

The recently developed Cgo field effect transistor (FET) devices (Schon et al. 2000a, b) claimed 
metallic and superconducting states for both electron and hole doping in the interface Cgo layer, 
the hole-doped system showing generally higher Tc than the electron-doped system. That could 
be related to a larger electron-phonon (e-ph) coupling of the HOMO-derived band than for the 
LUMO-derived tiu band (Manini et al. 2001). 

However, in an orbitally degenerate system like the one at hand, the electron-phonon coupling 
competes against intra-molecular exchange of Coulomb origin, responsible for Hund rules. In 
fact, Hund rules generally favor high spin for a degenerate molecular state, whereas coupling to 
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intra-molecular vibrations leads to a Jahn- Teller (JT) splitting of the degeneracy which favors low 
spin. Furthermore, since doped fullerenes are narrow-band molecular conductors, knowledge of the 
local Coulomb repulsion, usually parametrized by the so-called Hubbard U, is important in order to 
establish whether C^^ conductors are weakly or strongly correlated electron systems. In the former 
case, a conventional Eliashberg-type approach should be adequate to explain superconductivity. 
In the latter, a new theoretical framework is most likely needed (Fabrizio et al. 1997, Capone et 
al. 2000, Capone et al. 2002). 

All these considerations stress the importance of a realistic estimate of the Coulomb interaction 

terms (Hubbard U and the Himd multiplet terms) for Cg(^. In the past, electronic-structure-based 
calculations of these parameters have been made for the negative ions Cg,^ ' where the electrons 
are added to the t\u LUMO. In that case the structure of the Coulomb Hamiltonian is formally the 
same as that for an atomic p level in spherical symmetry, and as such entirely determined by two 
parameters only: the configuration-averaged U and Hund-rule intra-molecular exchange J (Martin 
and Ritchie 1993, Han and Gunnarsson 2000). For positive ions Cg^^, where n holes are added to 
the hu HOMO the only estimate available for the Coulomb parameters is a recent empirical one 
(Nikolaev and Michel 2002). The task of an electronic structure-based first principles calculation 
of these parameters will be the main purpose of the present paper. 

In a fivcfold-dcgcncratc hu orbital, as wc shall detail below, icosahcdral symmetry determines 
these Coulomb couplings in terms of five independent parameters: a configuration-averaged U, 
plus four intra-molecular exchange terms. All the low-energy electronic degrees of freedom of a 
solid-state system of positively (or negatively) charged Ceo molecules can be well described by a 
model Hamiltonian including the five hu hole bands (or three t\u electron bands) only, all other 
orbitals acting, as usual, as a mere source of renormalization of the 'bare' parameter values (Han 
and Gunnarsson 2000). 

In order to calculate the five independent intra-molecular electron-elcictron (e-c) Coulomb pa- 
rameters, wc use standard density-functional electronic-structure calculations in the local (spin) 
density approximation [L(S)DA], imposing as a constraint different values of the electronic occu- 
pation number in the individual Kohn-Sham (KS) orbitals. Wc carry out several frozen-geometry 
single-molecule constrained-LDA calculations of the total energy for a variety of charge and spin 
states of undistorted icosahedral Gqq. By comparing these energies with the corresponding an- 
alytic expressions for the model Hamiltonian, which is expressed in terms of the five unknown 
Coulomb parameters, we finally determine all of them. 

As the calculations are carried out for an isolated molecule, the computed Coulomb parameters 

are effective vahies, which contain the screening due to the polarizability of the filled molecular 
orbitals in the molecule, but contain neither the screening due to the other molecules nor that 
of the conduction electrons in the solid. As a reliability check, we also recompute with the same 
method the e-c U and J parameters for the LUMO band. The results are found in good agreement 
with previous estimates (Martin and Ritchie 1993, Han and Gunnarsson 2000, Antropov et al. 
1992), which further confirms the viability of our method. 

With the Coulomb parameters in hand we can then compute the multiplet spectrum for any 
given molecular occupancy, n =1,...5. This spectrum strictly applies only to ideal rigid Ceo ions, 
and is not of direct experimental relevance, because it leaves JT distortion effects out. The latter 
can be crudely estimated using the hole- vibration couplings previously calculated in Cgo (Manini et 
al. 2001) either at second-order in perturbation theory, corresponding to a full neglect of retardation 
effects, the so-called 'anti-adiabatic' approximation, or in the opposite 'adiabatic' limit. In the anti- 
adiabatic approximation, the effective e-e interaction is simply the superposition of the Coulomb 
repulsion and the phonon-mediated attraction. The total net result as far as U is concerned is still 
repulsive, the large Coulomb term only marginally corrected by molecular distortions. All the other 
intra-molecular exchange terms are instead heavily reduced. However, while in the case of Cg,^ 
that leads to an effective sign reversal from Hund to 'anti-Hund', (and from repulsive to attractive 
for an electron pair in the singlet channel) the balance is much less definite for Cg^, where the 
overall sign remains positive for n=2 and is uncertain for higher n values. The efficiency of the JT 
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effect in reversing Hund-rules couplings for CgQ*' is even weaker when the adiabatic approximation 
is considered instead of the anti-adiabatic one. In the adiabatic approximation, where ionic motion 
is classical, the molecular ground state of Cq^ turns out to be always high spin for all n values, in 
contrast to CgJ~ where it is always low spin. 

This paper is organized as follows: Sec. |^ introduces the model, determining the minimal number 
of independent parameters consistent with icosahedral symmetry. The constrained-LDA calcula- 
tion and its results are described in Sec. ^ The multiplet spectra resulting from the computed 
couplings are shown in Sec. |[ The results are discussed in Sec. ||, and some lengthy formulae are 
collected in an Appendix. 



2 The model Hamiltonian 

Our final target is to address the low-energy properties of (a lattice of) charged Cgo molecules. To 
this end we construct a model Hamiltonian to describe the physics of either the /i„ HOMO (holes) 
or tiu LUMO (electrons) bands. The role of the other orbitals is to act as a renormalization of 
the effective parameters for the band at the Fermi energy. In this paper we concentrate on the 
determination of single-molecule properties, and defer to a future work the calculation of the bands 
in the solid. The model Hamiltonian for a single molecule reads 



where 



H — Ho + Hvih + -ffc-vib + Hc-e (1) 



Ho = 4mCo-m (2) 



am 



H„, = ^^(/^^^+Q2^^) (3) 

iAfl 

-"e-vib — 2 ^-^ ^mm' V«A// ^am^f^rn' y"^) 

r iA a7nm' fi 

^e-e = ^YYw^^a'im,m';n,n')cl^cl,^,Ca'n'Can- (5) 

(7, a' mm' 

are respectively the single-particle Hamiltonian, the vibron contribution (in the harmonic approxi- 
mation), the electron- vibron coupling (in the linear JT approximation) (Manini et al. 2001, Manini 
and De Los Rios 2000), and the mutual Coulomb repulsion between the electrons. The c'l ,^ de- 
note the creation operators of either a hole in the HOMO or an electron in the LUMO, described 
by the single-particle wave function ipmai''^)- f indicates the spin projection, m and n label the 
component within the degenerate electronic HOMO/LUMO multiplet, and i counts the phonon 
modes of symmetry A {2 Ag, 6 Gg and 8 H^ modes). are Clebsch-Gordan coefficients of the 

icosahedron group, for coupling two hu (holes) or tiu (electrons) states to phonons of symmetry 
A. r is a multiplicity label, relevant for Hg modes only (Manini et al. 2001, Butler 1981). QiA/i 
and PiAp, are the molecular phonon coordinates and conjugate momenta. Spin-orbit is exceedingly 
small (Tosatti et al. 1996), and it is therefore neglected. 

The Coulomb matrix elements are defined by: 

where Ua,cr' (r, r') is an effective Coulomb repulsion, screened by the other electrons of the molecule. 

One way to estimate these matrix elements is to evaluate the Coulomb integrals (||) directly 
for the simple kernel UCT,CT'(r,r') = gg/(47reo|r — r'|) and given molecular orbitals (Nikolaev and 
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Michel 2002). This approach neglects completely the screening due to the other electrons on 
the same molecule. Here we choose a rather different approach: namely, we parametrize the 
interaction Hamiltonian (|^) in the most general way allowed by the molecular symmetry, and 
then determine the parameters by fitting to ab initio electronic structure calculations. As these 
calculations allow for the full polarization response of the total charge density (except for core 
levels, whose polarizability is negligible by comparison) the screening effect of all molecular valence 
electrons is accounted for in the final parameters. 

The symmetry of the Coulomb interaction plus the molecular symmetry of the problem allow 
us to express all of the Coulomb integrals in (|) as functions of a small set of physical parameters. 
In the following, we obtain the minimal set of Coulomb parameters that determine the interaction 
Hamiltonian (||), as required by the symmetry of the molecule and the symmetry label of the 
orbitals under consideration. 

As the Hamiltonian is time-reversal invariant, the orbitals can be chosen real with no loss of 
generality. Furthermore we take the orbitals, as well as the interaction, to be spin-independent 
(thus neglecting spin-dependent screening effects which might be possible in magnetic states), so 
that: 

Wcr,a'im,m';n,n') ^ w{m,m' -,71,71') = /d^r /d^r' 93„(r) (r') u(r, r') .^„(r) .^„/ (r') . (7) 



With the above assumptions one finds immediately: 

w{m, m ; n, n) — w{n, m! ; m, n) — w{m, n ; n, m) = win, n ; m, m') . (8) 

The effective screened interaction shows the full molecular symmetry, i.e. 

u(Rr,Rr') = M(r,r') (9) 

for all the symmetry operations R of the icosahedral group Ih- In order to make use of this symme- 
try, we decompose the product wave functions into irreducible representations of the icosahedral 
group: 

^„(r)^„(r)= ^ C;i^$;;^(r) , (10) 

using again the Clebsch-Gordan coefficients C,'J^j^ to couple two /i„ (holes) or tiu (electrons) tensors 
to an irreducible tensor of symmetry A. The label A runs in principle on all the irreducible 
representations (Ag, Tig, T2g, Gg, Hg, Au, Tiu, T2„, Gu, Hu) of the icosahedral group Ih- The 
multiplicity label r distinguishes between multiple occurrences of the same representation in the 
coupling (|lO|): it is the standard extra label for groups, such as It, which are not simply reducible 
(Butler 1981). Due to the symmetry relation (^), only the symmetric couplings occur. In particular, 
for holes in the HOMO, from the decomposition of /i„ ® the only nonzero contributions come 
from A = Ag, Gg, h'[=^\ and For electr ons in the tiu LUMO, we have K — Ag and Hg 

only. In terms of this symmetry recoupling, we rewrite the Coulomb matrix elements as: 

w{m,m';n,n') = ^ Clt' cUf"' U\ /dV $;;^(r) u(r, r') (r') . (11) 

r,A,fi 
r',A',fi' 

This equation shows the decomposition of the w{m,m' ;n,n') interaction matrix into the sum of 
products of geometric factors (Clebsch-Gordan coefhcients), times a relatively restricted number 
of coupled matrix elements. We can now exploit the symmetry of the Coulomb interaction (^ to 
further simplify the remaining integrals. To this end, we apply a generic group operation R to the 
integration variables. The explicit transformation properties of the coupled wave functions <f>J^^(r) 
allows to introduce the group representation matrices (R), while the effective interaction 



4 



remains invariant. Next, we apply the grand orthogonality theorem of representation theory, to 
rewrite the interaction as follows: 



dV y"dV$;^(r)u(r,r')«'P'^'(r') 

= Jd^r y"d3r'$[,'^(R-ir)M(R-ir,R-ir')$;;>'(R"^r') 

= 5]r;l^^(R)r;i;^,(R) Jd'r |dV'$;t(r)«(r,r')<f' (r') 



= Sa,A'S,,,>^^J2 /d'^ y'dV$;A(r)^(r,r')<'^(r'). (12) 

This shows that the integrals in Eq. ( |ll|) are diagonal in the representation label (but not in 
the multiplicity label r). This relation determines explicitly the most general expression for the 
Coulomb matrix elements in a shell of ti„ or icosahedral label, in terms of a minimal set of 
independent parameters F^'^ [defined in Eq. (|T^)]: 

w{m, m'; n, n') = ^^''''■^ f E C^") " (^3) 

r,r',A \ / 

In this paper we label states within the degenerate representation using the C5 quantum number m 
from the //i D D C5 group chain (Butler 1981). Note however that the purely geometric decom- 
position of the Coulomb integrals (^3|) holds for any choice of the group chain, and correspondingly 
of the Ih Clebsch-Gordan coefficients. 

In the case of electron doping in the tiu orbital, no multiplicity r labels appears, and thus, 
according to Eq. (|l3|), the Coulomb Hamiltonian is expressed as the sum of two terms, whose 
strength is governed by the two parameters F'^f and F^n. These parameters are related to the 
fc = and k — 2 Slater-Condon integrals i^*^*^-* for p electrons in spherical symmetry (Cowan 1981). 
For hole doping in the hu HOMO, we need five parameters 

Fi = F-^s, F2 = F'^\ F3 ^ F^^i'-^s, Fi = F'^'^'"^, F5 = F^-"^^"^ (14) 

to determine completely the Coulomb matrix elements Q In terms of spherical symmetry for a d 
atomic state, again corresponds to the totally symmetric F'^f parameter, while the F^^^ and 
F^'^'> spherical parameters are replaced by the four icosahedral F2^5. 

Rather than the F^ parameters, for ti„ electrons it is more common to use the parameters 
U = F^f /3 - F^' /3, and J = F^f /2. With this definition of [/ the multiplet -averaged energy 
has the simple dependence on the total number of electrons n: 

E-^%n) = Tr(i7|„) = en + C/I^il^ , (15) 

where H\n is the Hamiltonian restricted to the n-electrons states. The J parameter controls the 
multiplet exchange splittings, so that the center of mass of the multiplets at fixed n and total spin 

^The Coulomb Hamiltonian for icosahedral ?t-states was expressed in terms of five parameters also by Oliva 1997, 
Plakhutin and Carbo-Dorca 2000. 

It should be noted that U differs from the usual definition of the Hubbard U , involving the lowest multiplet 
in each n-configuration: f/™"^ = i?™"^(n -|- 1) -|- _E"""(n — 1) — 2i?™'"(n). This second definition is unconvenient, 
especially in the case, since it depends wildly on n. 
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S is given by |^ 



For the hu holes, the center of mass of the muhiplets with n holes is located at energy 



5 45 9 9 
This leads to the definition of an average Coulomb repulsion 



U 



Fi 4 F2 -F3 -F4 
"5 45 9 9" 



(16) 



(17) 



(18) 



We can define a spin-splitting parameter J also for the holes, by considering the center of mass of 
the multiplets at fixed spin S, £;^^<=(n, S). We find that the Coulomb Hamiltonian (|) is consistent 
with 

5 1 



with 



(19) 



(20) 



In what follows we take as a convenient set of independent Coulomb parameters: U, F2, F3, F4, 
and F5. 

Finally, with the decomposition ( [l3| ) in hand, it is convenient to re-organize the interaction 
Hamiltonian (||) in terms of number-conserving symmetry-adapted fermion operators: 



rr'A \ ^ 



-At 



where we defined the operators: 



w 



(21) 



(22) 



a mn 



and the constant A, which is [^U + | j) for holes and (^^U + 2J) for electrons. 



3 Determination of the Coulomb parameters 

After the explicit derivation of the form of a general icosahedral e-e interaction -ffc-o, we come now 
to the numerical calculation of the parameters fixing the interaction for Cgo ions. We compute 
these parameters by comparing the (analytical) expressions for the energies in the model Hamilto- 
nian (j^) with numerical results, obtained by first-principles density-functional theory (DFT) LDA 
calculations of the electronic structure of the Ceo molecule. 

As the Coulomb Hamiltonian governs the spectrum of multiplet excited states of the ionized 
configurations, in principle it would be straightforward to obtain the Coulomb parameters by 
fitting the excitation energies of (0) to multiplet energies obtained with some ab-initio method. 

In the tin case the eigenenergies can be written in the closed form 



5(5 + 1) + i-L(L + 1) + jn{n - 6) 



as function of n, 5 and the total 'angular momentum' L (recall that the tiu orbitals behave effectively as p-orbitals). 
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or to experimental multiplet spectra, if they were available and if the e-ph contribution could be 
separated out. However, to our knowledge no such experimental data are till now available. We thus 
choose to extract the Coulomb parameters from DFT calculations. Yet this is not straightforward, 
since in standard DFT the excitation energies of the KS system do not have a rigorous physical 
meaning; the KS states are only auxiliary quantities (Perdew 1985). Although excitation energies 
are accessible in DFT within the framework of time-dependent DFT (Petersilka et al. 1996), here 
we follow an alternative approach which is similar in spirit to the constrained-LDA (Dederichs 
et al. 1984, Gunnarsson et al. 1989) method to extract effective local Coulomb parameters. In a 
nutshell, constrained-LDA yields the ground-state energy of the system subject to some external 
constraint, such as a fixed magnetization or a fixed orbital occupancy. In practice, it is convenient 
to impose constraints such as to select states which are single Slater determinants, since they are 
described fairly accurately by standard LDA methods. By comparing their total energies with 
the expectation values of the model Hamiltonian with respect to the same states, it is possible 
to determine the interaction parameters, in the spirit of the ASCF (self-consistent field) scheme 
(Jones and Gunnarsson 1989). 

In order to describe the method, it is convenient to consider a simple example. Let us focus on 
the states n^; 0, 0; 0, 0) in the LUMO subspace where spin up and spin down electrons 
fill orbital (pi , the other two orbitals being empty, and define {n^^ , rij ) the corresponding total 
energies. One could determine the 'Hubbard Ui relating to this orbital as 



Ui = £;'°*(i, 1) + £;*°'(o, 0) - £;'°*(i, 0) - £;*°*(o, i) = £;'°*(i, i) + £;'°'(o, o) - 2e^°\i, o) . (23) 



A slight complication to this simple approach is brought in our problem by the orbital degeneracy. 
Suppose we indeed compute by LDA the energy of |1, 1; 0, 0; 0, 0). This orbital, where we place 
the two electrons, is initially degenerate to the other two LUMO components. However, in LDA, 
a Kohn-Sham (KS) filled orbital shifts immediately up in energy with respect to empty ones. 
Consequently, if we insist to fill one of the three, originally degenerate, KS orbitals, at the next 
iteration the electrons go naturally to occupy either of the two initially empty orbitals, with a 
large jump in charge density from one iteration to the next. This effect is due to the imperfect 
cancellation of the self interaction within LDA n, which is also related to the well-known LDA gap 



A possible remedy is to artificially introduce small gaps (of the order of 10 meV) in the otherwise 
degenerate HOMO and LUMO, by adding a tiny distortion of the icosahedral Cgo molecule along 
one of the JT active modes. Since the distortion is very small (each atom moving by less than 0.5 pm 
from its equilibrium position), the model remains essentially representative of fully symmetric 
fuUerene. However, in order to effectively control the charge (either electrons in the LUMO or 
holes in the HOMO) in the different orbitals, the splittings must overcome the self-repulsion. This 
fact suggests occupying an orbital by a small fraction of an electron so that the self repulsion is 
sufficiently small to leave this orbital in the same energy position dictated by the distortion field. 

To recover the actual value of the energy at integer charge, as required for the determination of 
U according to Eq. (p^), we make use of a known artifact of the LSDA: the ground-state energy 
of a system as a functional of the fractional occupation of a KS orbital interpolates smoothly the 
energies of integer multiples of the elementary charge. To the extent that the modifications of 
orbital (pi due to changes of its filling ni may be neglected, the total energy is a parabolic function 
of charge. In particular, in our specific example. 



where, according to Janak's theorem (Janak 1978), the linear coefficients b = b' equal the KS 
single-particle energy ei at ni — > 0. The quadratic coefficients c and c' are not generally the same, 

* Self-interaction is not the only possible origin of the convergence problem. The subtle problem of pure-state 
ii-representability leads to basically the same symptoms. See (Schipper et al. 1998). 



problem. 




(24) 
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since they are extracted from unpolarized and spin-polarized configurations, respectively. We have 
verified, in those configurations where the self interaction causes no convergence problem, that the 
extrapolation from calculations at ni < 0.2 is in very good quantitative agreement with direct 
total-energy calculations at integer rii's. Through Eqs. ( p^ ) and (|2^) we have: 

Ui = (^El°* + b-2 + y2^^+El°^~2 (^El°* + b' + |-f 2^ = (2c - c') . (25) 

We have therefore expressed the Hubbard Ui as a linear combination of the quadratic coefficients c 
and c' of the extrapolation parabolas. Note that the c's for both configurations involved are needed 
for the determination of C/i: it would be incorrect to identify Ui to, for example, the curvature c 
of the total LDA energy as a function of charge. 

The complete determination of the two and respectively of the five Coulomb parameters for ii„ 
electrons and for /i„ holes follows the same track as the simple determination of J7i outlined above. 
First, we select two sets of electronic configurations, one for electron and the other for hole doping, 
containing three (see Table |l|) and eight (see Table |^) elements, respectively. For each configuration 
we compute the total energy as a function of the (fractional) charge, for five values < rii < 0.2. 
The calculations for a given set are carried out with a fixed JT distortion within DFT-LSDA. As 
in previous calculations (Manini et al. 2001) we use ultrasoft pseudopotentials (Vanderbilt 1990) 
for C (Favot and Dal Corso 1999). The plane- waves basis set is cut off at 27 Ry (charge density 
cutoff = 162 Ry). The Ceo molecule is repeated periodically in a large simple-cubic supercell 
lattice of side a. To insure total charge neutrality, thus correcting for the G = Q divergence of the 
total energy, a compensating uniform background charge is added. The total energy is corrected 
for the leading power-law Coulomb interactions among supercells, by removing the Madelung 
term and the correction, with the method devised by Makov and Payne (1995). We extract 
the finite-a corrections by running several calculations with a ranging between 1.32 and 1.85 nm, 
as illustrated in Fig. |l| in a typical example. [It might have been marginally cheaper to use the 
modified Coulomb potential method (Jarvis et al. 1997) instead of the size scaling. That method 
however required a larger lattice parameter a, and thus more memory space]. Parabolas of the 
form ( p^ ) are fitted to the calculated a ^ oo energies. The resulting quadratic coefficients ci for 
electrons and holes are reported in Tables |l| and ^ respectively. 

In the light of Janak's theorem (Janak 1978), stating that the single-particle KS levels ei{ni) — 
dE^°^{ni)/dni, the Ci coefficients, besides representing second derivatives of the total energy w.r.t. 
charge, can alternatively be seen as first derivatives of the single-particle levels w.r.t. charge, as 
follows: 

e,(n,) = e,(0) + c, n, + 0{n^) . (26) 

From the a-scaling of the total energy (Makov and Payne 1995), we derive the a-scaling of the 
single-particle KS levels e^, which allows us to compute these quantities for the isolated molecular 
ion (a oo). Equation ( |2^ ) (neglecting 0(nf) corrections) provides a second method to derive 
the Ci coefficients, which are the basic ingredients in the calculation of the Coulomb parameters. 
As apparent in Tables |^ and ||, the coefficients Ci obtained from the total energy and from the 
single-particle levels are essentially in accord. However, the values from the single-particle levels 
are numerically more stable since, contrary to the total-energy method, they do not involve small 
differences of large numbers. In the following we shall use the c^'s from single-particle energies for 
the determination of the Coulomb parameters. 

The calculation of the e-e parameters is then realized by equating the various DFT extrapolated 
energies to the expectation values of (|l|) with respect to the same electronic configurations. Given 
the arbitrariness in the reference energy e for the model Hamiltonian (|l|) , we define for convenience 

e^U/2 + bi- fi . (27) 

where ^ is a free parameter allowing for the possibility of a chemical potential shift with respect 
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to the DFT calculation |. 

In Table | we collect the analytic expression of the energies of the three states considered 
for tiu electrons. Equating the terms on the third and fourth column of Table we have 3 
equations to determine the 3 unknown quantities U, J and fj,: we obtain the physical parameters 
simply by inversion of the linear dependency, and by replacing the values of Ci in Table |^. We 

have tabulated the combination 2n~^ (^{i\Hc^c\i) + Uni/2^ instead of simply {i\Hc-c\i), as each 

equation involves quantities of the same order of magnitude. By eliminating /i [in analogy to 



the one-state example of Eq. (25)], we find for the Coulomb parameters of the negative Ceo ions: 
U = fcin^n^Ti) - 5^10,1,0) = 3069 meV and J = |c|||^||^||) - |c||,|,|) + ^C|o,t^o) = 32 meV. These 
values, summarized in Table ^, are in the same range as previous estimates (Martin and Ritchie 
1993, Han and Gunnarsson 2000, Antropov et al. 1992). 

To produce a reliable estimate of the six unknown quantities (the five e-e parameters plus ^), 
we consider eight different /i„ hole states, for whose energies we collect the analytic expressions in 
Table |^. Therefore we have eight equations in six unknowns: we obtain the best estimate of the 
physical parameters by adjusting them to minimize the sum of the squared difference between the 
energies in the third and the fourth column of Table ||. For this overdetermined system of equations, 
the combination 2n~^ (^{i\Hc-c\'i) + U shows its advantage, that all LSDA calculations weight 

the same in the fit. This fitting procedure yields the values of the Coulomb parameters for CgQ*" 
collected in Table |[ The standard deviation of the fit (1.4 meV) gives an estimate of the numerical 
accuracy of the Ci parameters. By standard error propagation, we obtain the estimate of the 
errorbar on the individual e-e parameters reported in Table |^. 

We can now comment on our obtained results. We observe first of all that the only large 
parameter isU. It takes essentially the same value in the LUMO and the HOMO: this value of about 
3 eV governs the multiplet-averaged hole-hole repulsion, and is also compatible with experimental 
estimates (Antropov et al. 1992) for isolated molecular ions. In the solid, the screening of the 
local Coulomb parameters due to the polarizability of all the surrounding Cgo molecules could be 
approximately accounted for in a Clausius-Mossotti scheme (Antropov et al. 1992), and, for CgQ~, 
it may reduce U by roughly a factor 0.5 (Lof et al. 1992). The polarizability screening in the 
solid is expected to affect F2 ^ F5 much less than U. Note however that the actual Hubbard U, 
based on differences of ground-state energies, acquires an 71-dependent contribution of F2 F^. 
The appropriate [/™'" are collected in Table |^. The extra intra-molecular contribution is especially 
large at half filling. 

The relative smallness of the intra-molecular parameters compared to U is traced to the very 
close values of the quadratic coefficients Ci (listed in Table ||) for all different configurations 
In turn, this indicates that, contrary to the strongly localized orbitals of atomic physics, in Cgo it 
does not matter much the relative spin and orbital placement of two electrons in the HOMO or 
LUMO: they would always feel more or less the same repulsion of roughly 3 eV. The Hund rules are 
therefore rather weak in Cgo ions, because the degenerate orbitals are spread over a carbon shell 
of 7 Adiameter, rather than concentrated around a single nucleus. The largest parameter is F3 
corresponding to the r = 1 Hg symmetry. The F^ parameter is the smallest, effectively compatible 
within error bars with a zero value. 

The computed intra-molecular exchange J is almost twice as large ^ in the HOMO than in the 
LUMO (Tables ^ and Consequently, the splittings (in the order of hundreds of meV) of the 
n-hole Coulomb multiplets are larger for holes than for electrons. We come next to the detailed 
study of this multiplet spectrum. 



^Notice that /i is not exactly a chemical potential shift, since the linear term within our constrained DFT-LSDA 
has a component proportional to the Jahn- Teller splitting introduced to stabilize each set of electronic configurations. 

^ The definition of J for the ti^ orbital contains a factor 5/4 due to historical reasons, as apparent from the 
comparison of Eq. (^) and (p^. If this factor is accounted for, the actual ratio between the first-Hund-rule terms 
in the HOMO and in the LUMO is about 1.5 
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4 Multiplet energies 



The computed values of the couphng parameters can be used to calculate the multiplet spectrum 
of Hq-c- We concentrate here on the states of n holes in the fivefold degenerate /i„ HOMO. In 
order to diagonalize the Hamiltonian matrix, we wish to take full advantage of symmetry. For each 
charge n, (Ih D -D5) orbital label (A,//), multiplicity r, total spin S, and spin projection M, we 
first construct a set of symmetry-adapted states by iteratively coupling the one-hole hu orbitals to 
all the {n — l)-holcs states as follows: 

\n-l,T',A',fi';S',M') x |1, m"; ^, M") (28) 

where C^^^,.fi^ ^// and C'|,*^^,.^^2 m" icosahedral and spherical Clebsch-Gordan coefficients 

taking care of the orbital and spin recoupling respectively, and t' is the parentage of the (n — 1)- 
particle state. The resulting states \n, (r'. A', S'), (A, r),ii; S, M) (for all possible r'. A', S' and r) 
are then orthonormalized to form the set of n-hole basis states |n, r. A, ^; S', Af), with r counting 
their parentage. 

In this symmetry- adapted basis the Hamiltonian is diagonal with respect to the labels n. A, /i, S 
and M and its eigenvalues are independent of /i and M . This block-diagonal form allows in many 
cases to compute the analytical expressions for the multiplet energies -Emuit, i-e. the eigenvalues 
of Hc-c — Uh{n — l)/2, that are collected in Table ||. For 3 < n < 7, the calculation of the 
multiplet energies involves the diagonalization of block matrices of size 3 up to 7, where analytical 
methods are unpractical. For these cases we report the Hamiltonian submatrices for those states 
in Appendix |^. The analytical 3-holes spectrum shows the degeneracy of a Ti^ and a doublet 
state, which has been observed and explained in previous work (Oliva 1997, Lo and Judd 1999, 
Plakhutin and Carbo-Dorca 2000). 

The spectrum obtained by substituting the computed parameters of Table ^ into the expressions 
of Table ^ is collected in Table ||. For all values of n, we of course verify Hund-rule behavior, i.e. the 
high-spin state has the lowest energy. However the first Hund rule leads to comparable splittings 
to second Hund rule, so that states of different spin are energetically inter-mixed, for n = 4 and 5. 
The multiplet structure is similar to what was reported in Ref. (Nikolaev and Michel 2002), with 
a few differences in the detailed ordering of closely spaced levels. The total spread in the DFT 
results of Table ^, however, are a factor three smaller than those of that rigid-orbital unscreened 
calculation (Nikolaev and Michel 2002). 

The computed spectrum of Table ^ represents that of ideal rigid icosahedral Cgo- The coupling 
of electronic state to the intramolecular vibrations generally leads to JT distortions, involving 
energy scales that compete with the e-e repulsion and generally favor low-spin states, a sort of 
anti-Hund rule. The interplay of the Coulomb and e-ph terms originates a complex pattern of 
vibronic multiplet states that was studied in detail in the simpler case of the negative Cgo ions 
(O'Brien 1996). The HOMO system at hand is much more intricate, due to the interplay of several 
parameters. Here we shall address this problem at a more approximate level of accuracy. 

First we observe that, in the limit where the typical phonon energies hiVi are much larger than the 
e-ph energy gains ~ gf'hijJi, the e-ph Hamiltonian treated at second order in gi takes the form of the 
first term in the right hand side of Eq. ( [2l| ) [anti-adiabatic or weak-coupling limit). The strengths 
of the effective e-e interaction parameters are given in terms of the dimensionless couplings 5]"^ 
and frequency cj^a by F^'^ — — J2i 9iA9iA^^i^/^ Q Table || hsts the parameters Fi generated 

According to (Manini et al. 2001), the coefficients for tfie coupUng of Hg with Hg to Ag and to Gg are normahzed 
so that ^jj^ ^C'nm''^ = 5 and ^^^^ {^nm'^^ = 5/4- Here we prefer to apply the standard normalization to 

unity, thus we include the 5 and 5/4 factors into the coupling parameters. Within this convention, a factor 6 
must be included in for the coupling of t^u electrons to Hg modes. 
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by the e-ph, in the notation of Eqs. ( p| ) (|T8|) and (|2^), compared to the Fi Coulomb parameters. 
At this level of approximation, the e-e and e-ph terms are expressed as smiis of formally identical 
terms, differing only for the value of the parameter multiplying each term. Thus, it is natural to 
combine the two contributions into a total effective two-body Hamiltonian which is formally 

identical to the definition of Eqs. (H) and ( p^ but based on total effective e-e parameters given by 
the algebraic sum of the Coulomb and e-ph parameters. These total effective e-e parameters are 
listed in the last column of Table 

The contributions to U of the Gg and Hg phonons [Eq. (p^)] is larger than the Ag term, thus 
giving a small repulsive phononic multiplet-averaged U. However, the huge Fi Coulomb 'monopole' 
parameter is barely affected by the Ag phonons-originated term. Accordingly, the total multiplet- 
averaged interaction J7*°* is left basically unaffected by the phonons contribution, for both the 
HOMO and the LUMO. 

On the contrary, the intra-molecular Hund terms Fi are of values comparable to the phononic 
counterpart, thus leading to a strong cancellation. In particular, for the holes the largest repulsive 
term F^ is reversed by the even larger coupling F^ to the Hg^~^^ phonon modes. The largest 
total effective term is F2: it remains positive, due to the modest Gg-phonons-mediated attraction. 
The corresponding multiplets spectrum, is reported in Table and drawn for n = 2 and n = 3 
in Fig. 1^. For n — 2 holes the ground state remains a triplet, with a very small gap to the 
lowest singlet, while low-spin ground states prevail for n > 3. For electrons, the total effective 
J = —25 meV indicates that low-spin anti-Hund states are to be expected for CgQ~, as is indeed 
observed experimentally in many electron-doped Ceo compounds (Kicfl et al. 1992, Zimmcr et al. 
1995, Lukyanchuk et al. 1995, Prassides et al. 1999). 

The anti-adiabatic approximation, although tending to overestimate the e-ph energies is ad- 
vantageous in allowing to map the e-ph Hamiltonian onto an effective attractive e-e term: this 
mapping does not apply any longer when the coupling energies are not taken as much smaller than 
the harmonic phonon energies. However, the relatively large values of the realistic e-ph coupling 
in both positive and negative Cgo ions make the weak-coupling approximation not truly justified. 
Indeed the e-ph energetics based on this approximation are grossly overestimated if applied to 
intermediate/strong couplings. In practice, the JT energy gains and gaps in units of gffiuji become 
significantly smaller as the coupling changes from weak {gi <C 1) to strong {gi ^ 1). In particular, 
for tiu electrons interacting only with Hg modes (no e-e terms), the energy lowering divided by 
gfhoji drops by 60 % from weak to strong coupling (Manini et al. 1994). For hu holes, for n = 2 and 
a single Hg~'^^ mode (O'Brien 1972, Manini and Tosatti 1998), we see in Fig. ||that the spin gap 
in units of g^fiu reduces by 17 % only, going from the weak to the strong-coupling limit. However, 
when all Ceo modes are included, this reduction is as large as ^50 % due to contributions of the 
hI^~'^'^ and Gg modes. 

As the actual (anti-Hund) e-ph coupling should have weaker effects than those estimated in the 
anti-adiabatic approximation, the question of what is the symmetry of the ground state of the Cg^*^ 
ions remains open. The case of C^J marks an exception, since Coulomb couplings prevail already 
at the anti-adiabatic level (Fig. ||). Thus the prediction of an S* = 1 magnetic ground state for 
n=2 holes in Ceo seems fairly robust, at least within LDA accuracy. In order to settle this problem 
for the other cases, we study the e-ph coupling in the opposite, adiabatic limit, which becomes 
exact in the limit of strong e-ph coupling, and which proved quantitatively more realistic for Cg^*^ 
ions (Fig. ||). At the adiabatic level, the phonons are treated classically, with the electrons/holes 
contributing through (||) to the total adiabatic potential acting on the phonon coordinates QiA/i 
which are treated as classical variables. The additional ingredient we include here, and which was 
not included in the previous adiabatic calculation (Manini et al. 2001) is the e-e coupling. In the 
JT-distorted configuration, the icosahedral symmetry is broken, therefore all symmetry states are 
mixed. Only n, S and Sz are conserved. For example, the coupling to the distortion of n = 3 
holes, S = 3/2, M = 3/2 mixes the 10 states of Ti„, T2„ and T2U symmetry listed in Table ^. 
Assuming that the Coulomb parameters Fi are unchanged upon distortion, for each n, S and M 
we determine the optimal distortion, by full minimization of the lowest adiabatic potential sheet in 
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the space of all the plionons coordinates. In Table |ll|, we report the resulting lowest-state energy 
in each spin sector, based on the e-e and e-ph couplings of Cg^'^ ions. As already stated, we find 
a difference between electrons and holes: the ground state is low spin for Cg^ electrons, while it 
is always high spin for Cgg*". The case of Cg^ has almost degenerate 5 = 1/2 and S = 3/2 states, 
the former probably prevailing once non-adiabatic corrections are accounted for (Capone et al. 
2001). For the positive ions instead, the adiabatic result overturns the anti-adiabatic prediction 
of low-spin ground state for n > 3 holes. In these cases, the adiabatic spin gaps are fairly large, 
and are likely to survive when zero-point quantum corrections to the adiabatic approximations are 
added. For Cqq the gap to the lowest singlet state is rather small, but, as noted above, here the 
ground state is a spin triplet even in the anti-adiabatic approximation: it is likely to remain high 
spin also within an exact treatment of the phonons. 

The outcome of the adiabatic calculation is that positive Cgo ions favor high-spin ground states, 
while in negative ions e-ph coupling prevails and low-spin ground states are likely. However, for 
Cg^'^ ions the balance of e-e and e-ph is rather delicate, therefore the problem of the spin symmetry 
of the ground state ions is far from trivial, and remains basically open. Indeed, in some chemical 
environments high-spin states are observed to prevail in negative fuUerene ions (Schilder et al. 
1994, Arovas and Auerbach 1995), and this indicates that the lowest niultiplets of different spin 
type are almost degenerate. To get a more conclusive answer on this point for ions of both signs, 
it would be crucial to carry out a full diagonalization of (|l|) including all the phonon modes and 
Coulomb terms on the same ground, on the line of O'Brien (1996): we plan to carry out such 
calculation in a future work. 

5 Discussion and Conclusions 

The Coulomb couplings of holes in Cgo obtained in this paper are based on rigid icosahedral ge- 
ometry calculations. However, clearly in each different-charge state, the Cg^*^ molecular ion relaxes 
to different equilibrium positions, according to the interplay of e-ph coupling {Hvih + ^^o-vib) with 
intra-molecular Coulomb exchange (ffc-c)- In principle one could compute the Coulomb param- 
eters, allowing simultaneously for geometry relaxation. The disadvantage of such a calculation is 
the difficulty of disentangling the e-ph and e-e contributions. A second difficulty of principle is that 
the ion, in an electronically degenerate state, distorts to several equivalent static JT minima of 
less than icosahedral symmetry (Manini and De Los Rios 2000). These local minima are connected 
by tunneling matrix elements which mix them to suitable dynamical combinations of the different 
distortions, thus restoring the original icosahedral symmetry: such non-adiabatic situation would 
be outside the range of applicability of current standard first-principles computational methods, 
usually based on the Born-Oppenheimer separation of the ionic and electronic motions. Moreover, 
the lack of exact cancellation of self-interaction in the LDA makes even a practical attempt at a 
static, adiabatic calculation for JT-distorted ions impossible at this stage. These are the reasons 
that suggested restricting this first Coulomb calculation to the rigidly undistorted geometry. 

A further limitation of the present calculation is the assumption that the Coulomb parameters 
are independent of the charge of the state. In principle, due to both orbital and geometrical 
relaxation, the effective Coulomb interaction (^) will depend on the instantaneous charge state of 
the fuUerene ion. However, this effect, a very important one in single-atom calculations, is expected 
to be fairly small in such a large molecule as Cgo . Thus our parameters represent an average over 
n. 

For Cgo^, the value of J was estimated in the 100 meV region by Hartree-Fock calculations 
(Chang et al. 1991), and by direct integration of the unscreened Coulomb kernel (Nikolaev and 
Michel 2002). Both these methods overestimate the Coulomb repulsion because of underestimation 
or complete neglect of screening. The LSDA, where we get J = 32 meV is on the other hand 
known to overestimate screening, and thus to underestimate the Coulomb parameters. Some value 
in between, such as J '-^ 50 meV, as suggested by Martin and Ritchie (1993) is probably a more 
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realistic estimate of J in Cqq . Coming to the C^q case, we can regard the e-e coupUngs of Table || 
as lower bounds, the actual repulsion being possibly a factor 1.5 or 2 larger. Indeed, the calculation 
of Nikolaev and Michel (2002) finds splittings about 3 times larger than those of Table ^ and those 
can reasonably be regarded as upper bounds. 

On the other hand, the competing e-ph interaction is also very likely underestimated by LDA, 
as was demonstrated in the case of C^q by the comparison of the calculated couplings to those 
extracted from fitting the photoemission spectrum (Gunnarsson et al. 1995), which suggested values 
roughly twice as large. In that case, the effective e-ph J, determined from the photoemission data 
is J = —127 meV compared to the LDA value of J = —57 meV. In conclusion, both the Coulomb 
repulsion and the phonon-mediatcd attraction calculated within LDA are likely to need a rescaling 
by a similar factor of order two. Thus the balance between the two opposing interactions remains 
delicate in both Cg^j" [as demonstrated by the presence of both high-spin and low-spin local ground 
states in different chemical environments (Brouet et al. 2001, Kiefl et al. 1992, Zimmer et al. 
1995, Lukyanchuk et al. 1995, Prassides et al. 1999, Schilder et al. 1994, Arovas and Auerbach 
1995)] and even more so in CgQ^, where however the high-spin states should be more favored. 
Moreover, especially in the hole doped case, we find multiplet splittings which are comparable 
to the theoretical bandwidth of solid-state fullerene (about 0.5 eV), indicating that Hund-rule 
intramolecular interactions are an important ingredient in Cgo ions. 

Since any treatment of superconductivity caused by the JT coupling must include the competing 
Hund-rule terms, our results surprisingly indicate that positively-doped Cgo could display a weaker 
tendency toward superconductivity than negatively-doped Cgo- Magnetic states could occur for 
any integer hole filling, more commonly than for integer electron fillings. Even if magnetism were 
to be removed owing to band effects, one should still expect Cqq to make better superconductors. 
This conclusion is unexpected in the light of recent data claiming a larger superconducting Tc in 
positively charged than in negatively charged Cgo FETs (Schon et al. 2000a, b). The reasons for 
this disagreement are presently unclear, and will require further theoretical and experimental work. 
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A Appendix 

We provide here the recipe to construct the matrices i?inuit['T', A, S*], whose eigenvalues i?muit give 
the {F2 . . . F5)-dependent contribution to the multiplet energies of n holes of global symmetry A 
and total spin S in Table ^. Each matrix H-^^\t[n^ K, S] is a linear combination of four numerical 
matrices Mi[n, A, S] (given below), with as coefficients the Coulomb parameters F2 . . .F^,: 

5 

il,„uit [n. A, 5] = ^ F, AU [n, A, S] . (29) 

i=2 

With specific choice of the parameters, it is possible to study the effect of one particular operator: 
for example, by taking only ^3 7^ 0, thus diagonalizing the Al^[n, A, S] matrices one may study 
analytically the multiplet spectrum associated to that operator. To get the spectra of Tables ^ and 
we plugged the parameters of Table || into Eq. (p9|), and proceeded to diagonalize i/muitl", A, S] 
numerically. 
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The matrices Mj[n, A, S] are as follows: 
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Af5[5,if„,l/2] 



4 

4,A 

3 

9 /Is 
V 21 

4^/¥ 



4 

V 21 


4 

1/21 


\/2 
3 


9 /lO 
V 21 


4 ^/X 

* V 21 

3 


V 3 \ 

3 


8 vS 


2 \/2 


V 2T 


/T7T 

vio 


19 \f2 


2 \/^ 


21 


7 


3 


21 


63 


9 


2 ^2 


8 ^5 


5 V 2T 


1 


31 ^5 


13 


7 


21 


6 


14 


126 


18 V7 




6 









1 


~3~ 


6 


6 \/3 


\/To 

21 


1 

14 


3 


3 


7 

18 


4^/f 
9 


19^2 


31 V5 


\/¥ 


7 


4^5 


31 


63 


126 


6 


18 


9 


18%/7 


2 


13 


1 




31 


V5 

9 / 


9 


18^7 


6^3 


9 


18^7 



We report here the explicit expressions for the quantities A1...A5 of Table originated by 
diagonalizations of 2 x 2 blocks: 



Ai = i [4 i^l- 12 F2F3 + 9 F| + 4 F2^4- 6 i^3^4 + F| + 9 



A 



A2 
A3 

4/5 



^ [4 F| + 12 F2F3 + 9 F| - 20 F2F4 - 30 F3F4 + 25 + 300 



1 



— [4 F| + 36 ^2^3 + 81 F| - 44 F2F4 - 198 F3F4 + 121 F| + 288 FI\ 



21 1/2 



— [964 F| - 324 ^2^3 + 81 FI - 1604 ^2^4 + 162 ^3^4 + 721 F, 



T984\/5F2F5 ± 108 %/5F3F5 ± 876 %/5^4i^5 + 1332 



1/2 



(30) 

(31) 
(32) 

(33) 
(34) 
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State 


c^ [eV] 


c, [cV] 


\i) 


(from 


(from a) 


1 Ti,n,Ti) 


3.06855 


3.06850 


1 T,T,T> 


3.04652 


3.04659 


|o,T,o) 


3.06581 


3.06650 



Table 1: Quadratic coefficients c, for the different configurations used in the Cg,^ calculations, 
obtained from the total-energy curvature and from the linear dependence of the single-particle KS 
energies. The Slater determinant states \i) are specified by the individual occupancy of the LUMO 
orbitals labeled by m = — 1, 0, 1. 



State 


Ci [eV] 


Ci [eV] 


\i) 


(from E^°^) 


(from Ej) 


|o,o,Ti,o,o) 


3.20874 


3.20251 


|o,n,o,n,o) 


3.12273 


3.12234 


1 Ti,o,o,o,Ti) 


3.12349 


3.12274 


1 Ti,Ti,Ti.Ti.Ti) 


3.10285 


3.10090 


1 T,T,T,T,T) 


3.08025 


3.07928 


|o,T,T,T,o) 


3.07822 


3.08319 


1 T,o,T,o,T) 


3.08654 


3.08415 


|o,o,T,o,o) 


3.15757 


3.15165 



Table 2: Quadratic coefficients Ci for the different configurations used in the Cqq calculations, 
obtained with the two methods. The \n-2,n-i,no,ni,n2) notation specifies the occupancies of 
the m = —2, —1, 0, 1, 2 HOMO orbitals in each Slater determinant i. 
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State \i) 


Hi 


( 


'^\H,_,\i) + Unj2 


) 






Model: U + F[\i] 


] from Eq. (||): 




1 n,n,n) 


6 


U 






1 T,T,T> 


3 


U-lJ 


C|T,T,T) ^■ 




|o,T,o) 


1 


u 


C|o,T,o) ^ 


- 2/i 



Table 3: Comparison of the energy expectation values in the model (|l|) and from the LSDA 
extrapolation (p3) for the electron states considered in the calculation. 



Parameter 




(/fe-vib, 2"^^ order) Total: (e-e)- 


f (e-vib) 




[meV] 


[mcV] 


[meV] 


U 


3069 


32 


3101 


J 


32 


-57 


-25 



Table 4: The Coulomb parameters for Cgp , as obtained from the LSDA calculations, the effective 
parameters obtained from second-order treatment of _ffc-vib [based on the couplings of (Manini et 
al. 2001)], and the sum of the two contributions. 



State \i) rii 2n^ ^ i^{i\Hc-c\i) + Unj2'^ 

Model: U + F[\i)] from Eq. (p4|): c, + 2^i/ni 



o,o,Ti,o,o) 


2 


U + 


2 77' 

45^2 


o,n,o,Ti,o) 


4 


u - 


15^2 


n,o,o,o,Ti) 


4 


u - 


1^^2 


Ti,Ti,Ti,Ti,Ti) 


10 


u 




T,T,T,T,T) 


5 


u - 


45 2 


o,T,T,T,o) 


3 


u - 


27 2 


T,o,T,o,T> 


3 


u - 


27 2 


o,o,T,o,o) 


1 


u 





l^^3 + jgFi C|o,o,Ti:0,o> + /-i 

11 F _i_ 11 F V5 F ^ ^ 1 

240^3 + 48-^4 - —^5 C|o,Ti,0,Ti:0) + 

+ + C||j^^o,o,o,Ti> + 

C|n,Ti,Ti,Ti,n> + Ia* 



- iiFi + C|o,t,t,t,o) + ¥ 

F5 C||^o,T,o,T> + ¥ 

C|o,o,T:0,o> + 



Table 5: Comparison of the energy expectation values in the model (|l|) and from the LSDA 
extrapolation (EJ) for the hole states considered. 
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Parameter 




(i?e-vib, 2"d order) 


Total: (e-e)+(e-vib) 




[meV] 


[meV] 


[meV] 


Fi 


15646 ± 9 


-18 


15628 


F2 


105 ± 10 


-62 


42 


F3 


155 ± 4 


-173 


-18 


Fi 


47 ± 5 


-50 


-3 


F5 


± 3 


-14 


-14 




-27 ± 1 






U 


3097 ± 1 


27 


3124 


J 


60 ± 1 


-57 


3 



Table 6: The Coulomb parameters for Cgg", as obtained from the LSDA calculations described in 
the text. Two of the tabulated parameters (e.g. Fi and J) are linear combinations of the five others. 
The table also includes the effective (negative, anti-Hund) parameters obtained from second-order 
treatment of Hc-vih [based on the electron-phonon couplings of Manini et al. (2001)], and the sum 
of the two contributions. 



n 


^min 


spin 


1 


3038 




2 


3077 


76 


3 


3076 


99 


4 


3038 


132 


5 


3415 


202 



Table 7: The Coulomb Hubbard J/™^" = E"''"'{n +1) + E"''"'{n - 1) - 2E"''"'{n) and the Coulomb 
spin-gap Ag~;^ to the first low-spin state. 
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State 
symmetry 



spin Energy 

S -Emult 















71 = 


4 


















^9 1 


x3] 





Eigenvalues(i7,nuit [4, Ag,0]) 














Tig 







iF2 - j^Fs + iF4 - ^F5 


State 


spin 


Energ 


y 






T2g 







iF2 - j^Fs + \Fi + ^Fs 


symmetry 


S 










Gg 


[x4] 





Eigenvalues(77,nuit [4, Gg, 0]) 


n = 2 












Hg 


[x5] 





Eigenvalues(7?niuit [4, Hg,Q]) 


Ag 













Tig 


[x3] 


1 


Eigenvalues(iJmuit [4, Fig, 1]) 


Gg 





1 

— Fo - 

18^2 




25 7~i 
36-^4 




T2g 


r oi 

[x3] 


1 


Eigenvalues(iimuit [4, r2g, IJ) 


Hg 





1^2 + 


M-^3 + 


36-^4 ^ 


Ai 


Gg 


[x3] 


1 


Eigenvalues(iii„uit[4, Gg, IJ) 


Hg 





IF2 + 


36 3 


M^4 + 


Ai 


Hg 


[x3] 


1 


Eigenvalues(ffi„uit[4, ^^g, 1]) 


Tig 


1 


-1^2 


36^3 


+ ^^4 


+ ^F5 


Hg 




2 


2 TT" 5 77' 5 TP 

-^t2 - 6^-3 - 6-^4 


T2g 

Gg 


1 
1 


-1^2 
18^2 - 


- ^^3 

_ - 
36 3 


+ m^i 

23 p 
36-^4 




n = 
Au 


5 


1/2 


5 77' 1 13 77' 1 1 77' A 

^18-^2 + 36-^3 + 36^4 - A3 


n — 3 












A 

■^U 




1/2 


5 TP 1 13 TP 1 1 IT' 1 A 

"18^2+36^^3+ 36^4 + A3 


Tiu 
Tiu 


1/2 
1/2 


6 2 
i^^2 + 


12^ 3 
iFg - 


+ 15-^4 
iF4 


-4f. 


Tiu 
Tiu 




1/2 
1/2 


1 7-^ 1 1 3 7-^ 1 1 7-1 \/K 7-1 

— Ft -1- — Ft — — F/i — Ft; 

18-^2 -1- 3g-f3 3g-f4 2 

7-, ,17 7-1 , 1 7 7-n , '^^/^ 7-n A 

§F2 + i|F3 + i|F4 + 2^F5 - A4 


T2u 


1/2 


3^2 + 


i^3- 


iF4 




1^ 




1 /o 
1/2 


23 77' 1 17 77' 1 17 77' 1 3\/5 77' 1 A 

36-^2 + 72-^*3 + 72-^^4 + —i'b + A4 
18-^2 -1- 3g-f 3 3g-f4 -1-2-^5 


T2u 


1/2 


4^2 


12^ 3 


+ YiF4 


+ f F5 


T2U 




1/2 


Gu 


1/2 




+ — Ft 
^ 12^3 


+ j^F4 


-A2 


T2u 




1/2 


23 17 1 17 1 17 IT' 3v^ IP A 
3gi<2 + 72-^^3 + 72-^4 4~^^ ~ ^5 

23^ , lZi?„ + lift _ 3V5i^, + A, 
36 2T^72-'3T^72^4 4 ^5T^ "5 


Gu 


1/2 


-i^2 


+ — Ft 
^ 12^3 


+ S-^4 


+ A2 


T2u 




1/2 


Hu [x4] 


1/2 


Eigenvalues (Hmuit [3, 1/2]) 


Gu 


[x5] 


1/2 


Eigenvalues(iJmuit [5, G„, 1/2]) 


Tiu 


3/2 


-2^ 

3^2 


- ^^3 


~ T2Ti 


+ ^F5 


Ftu 


[x7] 


1/2 


Eigenvalues(iJinuit [5, Hu, 1/2]) 


T2u 


3/2 


_2^ 
3 2 


- ^F, 

12^ 3 


- T2T4 


Tiu 




3/2 


2 77' 5 77' 17 77' V5 77' 

9^2 - 36i^3 - 36-^*4 " —^5 


Gu 


3/2 


ii^, - 
6-^2 


^F, - 

12-^3 


11 p 

T2^4 




T2U 




3/2 


9^2 3g-f 3 3gJ^^4 2 ^5 














Gu 




3/2 


18 2 36-^3 36^4 














Gu 




3/2 


_ 11 _ A ^„ J_ 13 F, 
18-^2 3g-C3 -r 3gX'4 














Hu 




3/2 


_4^ _ 23^ _ lip _ A 
9 2 36-^3 3g-'4 "1 














Hu 




3/2 


_4^ 23^ 1 ^ 
9-^2 3g-C3 36-'4T^^l 














A-u 




5/2 


10 p 25 p 25 p 
"T ^2 - Tg^3 - 18^4 



Table 8: The Coulomb multiplets for Cgo+ as a function of the e-e parameters. The model 
Hamiltonian (|^) obeys particle-hole symmetry: therefore the multiplet energies for n > 5 holes 
equal those for (10 — n) holes. The non particle-hole symmetric contribution [en + Un{n — l)/2] 
is left out in this table. The H,nait[n, A, S] and A; quantities are defined in Appendix 
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n 


A 


s 


degeneracy 












1 





1 


Hu 


1/2 


10 





2 


T2g 


1 


9 


-59 




Tig 


1 


9 


-59 




Gg 


1 


12 


-11 




Gg 





4 


17 




Hg 





5 


46 




Hg 





5 


115 




^9 





1 


318 


3 


T2u 


3/2 


12 


-138 




Tiu 


3/2 


12 


-138 




Gu 


3/2 


16 


-91 






1/2 


6 


-39 




Tiu 


1/2 


6 


-38 




Hu 


1/2 


10 


-9 




Hu 


1/2 


10 


6 




Gu 


1/2 


8 


23 




T2u 


1/2 


6 


71 




Tiu 


1/2 


6 


71 




Gu 


1/2 


8 


96 




Hu 


1/2 


10 


99 




Hu 


1/2 


10 


247 


4 


Hg 


2 


25 


-238 




Gg 


1 


12 


-106 




Hg 


1 


15 


-94 




Gg 


1 


12 


-61 




T2g 


1 


9 


-59 




Tig 


1 


9 


-59 




Gg 





4 


-41 




Ag 





1 


-18 




Hg 


1 


15 


8 




^9 





5 


8 




T2g 


1 


9 


9 




Tig 


1 


9 


9 




Hg 


1 


15 


19 




Hg 





5 


34 




T2g 





3 


51 




Tig 





3 


51 




Ag 





1 


65 




Gg 





4 


113 




Gg 


1 


12 


120 




Hg 





5 


136 




T2g 


1 


9 


137 




Tig 


1 


9 


137 




Gg 





4 


146 




Hg 





5 


186 




Gg 





4 


256 




Hg 





5 


280 




Ag 





1 


494 



n A 


S 


degeneracy 


i^miiit fmeVl 




5/2 


6 


-397 


H^i 


3/2 


20 


-195 




3/2 


20 


-125 


Gil 


3/2 


16 


-97 




1/2 


10 


-82 




1/2 


8 


-69 




3/2 


16 


-68 


■^u 


1/2 


2 


-62 


T2u 


3/2 


12 


-21 


Tiu 


3/2 


12 


-20 


Gu 


1/2 


8 


-14 


Hu 


1/2 


10 


3 


Hu 


1/2 


10 


9 


T2u 


1/2 


6 


47 


Tiu 


1/2 


6 


47 


Tiu 


1/2 


6 


54 


T2u 


1/2 


6 


55 


Hu 


1/2 


10 


55 


Hu 


1/2 


10 


79 


Gu 


1/2 


8 


101 


A-u 


1/2 


2 


118 


Gu 


1/2 


8 


156 


Gu 


1/2 


8 


160 


T2u 


1/2 


6 


175 


Tiu 


1/2 


6 


175 


Hu 


1/2 


10 


185 


Hu 


1/2 


10 


334 



Table 9: The Coulomb multiplet energies for the Coulomb parameters obtained for Cg^. 
The states are sorted by increasing energy. The non particle-hole symmetric contribution 
[en + Un{n — l)/2] is not included in these numbers. 
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n 


A 


s 


degeneracy 












1 





1 


Hu 


1/2 


10 





2 


Tig 


1 


9 


-33 




Hg 





5 


-21 




T2g 


1 


9 


-1 




Gg 





4 


3 




Ag 





1 


15 




Gg 


1 


12 


21 




Hg 





5 


27 


3 


Hu 


1/2 


10 


-50 




Gu 


1/2 


8 


-45 




Tiu 


3/2 


12 


-36 




T2u 


1/2 


6 


-18 






3/2 


12 


-5 




Gu 


1/2 


8 


-4 




Hu 


1/2 


10 


-4 




T2u 


1/2 


6 


9 




Tiu 


1/2 


6 


9 




Tin 


1/2 


6 


14 




Gu 


3/2 


16 


17 




Hu 


1/2 


10 


49 




Hu 


1/2 


10 


58 


4 


Ag 





1 


-90 




Hg 





5 


-79 




Tig 


1 


9 


-54 




T2g 


1 


9 


-47 




Hg 





5 


-41 




Gg 





4 


-39 




G, 


1 


12 


-38 




Hg 


1 


15 


-31 




Gg 





4 


-18 




Hg 


1 


15 


-13 




Hg 


2 


25 


-11 




Tig 


1 


9 


-2 




T2g 


1 


9 


2 




Hg 





5 


3 




T2g 





3 


6 




Gg 


1 


12 


6 




Ag 





1 


10 




Gg 





4 


21 




Hg 





5 


23 




Hg 


1 


15 


31 




Tig 





3 


38 




T2g 


1 


9 


43 




Gg 


1 


12 


50 




Tig 


1 


9 


54 




Hg 





5 


63 




Gg 





4 


101 




Ag 





1 


116 



n A 


S 


dcgeneracv 


i^miiit fmeVl 


5 Tiu 


1/2 


6 


-84 


Hu 


1/2 


10 


-74 




1/2 


8 


-49 




1/2 


8 


-45 




1/2 


2 


-39 




1/2 


10 


-31 




3/2 


20 


-30 


Gu 


3/2 


16 


-24 




1/2 


6 


-19 




5/2 


6 


-18 


Hu 


1/2 


10 


-17 


Hu 


1/2 


10 


-8 


Gu 


3/2 


16 


-7 


T2u 


3/2 


12 


-2 


Au 


1/2 


2 


2 


T2u 


1/2 


6 


4 


Hu 


1/2 


10 


8 


T\u 


1/2 


6 


12 


Gu 


1/2 


8 


17 


Hu 


3/2 


20 


17 


Tiu 


3/2 


12 


29 


Gu 


1/2 


8 


33 


Hu 


1/2 


10 


41 


Gu 


1/2 


8 


48 


Tiu 


1/2 


6 


81 


T2u 


1/2 


6 


88 


Hu 


1/2 


10 


91 



Tabic 10: The Cgo+ multiplct energies for the effective Coulomb parameters, including the e-ph 

coupling in the anti-adiabatic (weak coupling) limit. The states are sorted by increasing energy. 

The non particle-hole symmetric contribution [en + U^°*"n{n — l)/2] is not included in these results. 
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n± 


S 


adiabatic 


anti-adiabatic 


2+ 





-129 


-21 




1 


-142 


-33 


3+ 


1/2 


-168 


-50 




3/2 


-222 


-36 


4+ 





-200 


-90 




1 


-211 


-54 




2 


-308 


-11 


5+ 


1/2 


-203 


-84 




3/2 


-256 


-30 




5/2 


-397 


-18 


2- 





-92 


-100 




1 


-71 


-25 


3- 


1/2 


-85 


-50 




3/2 


-97 


+75 



Table 11: The energy (in meV) of the lowest state for each n and S, including the e-e and e-ph 
contributions from _ffvib + ^^c-vib + ^^c-c (but excluding the [en + Un{n — l)/2] term), for hu holes 
(upper panel) and ii„ electrons (lower). First column: the phonons are treated in the adiabatic 
(strong-coupling) approximation, by full relaxation of the phonon modes to the optimal classical 
JT distortion for each n and S. Second column: the lowest energy for each n and S is reported 
from Table |l^. The anti-adiabatic energy lowerings are smaller than the adiabatic ones because of 
a larger cancellation of e-e and e-ph contributions. 
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Figure 1: The total LDA energy (diamonds) of the Ceo molecule as a function of the inverse 
lattice parameter a^^ for (uniformly spread) hole charge i ,t i ,u ,U ,U) — 0-2 (main graph) and 
""ITi Ti Ti.Tl.TD ~ ^ (inset). The solid lines are finite-size polynomial scalings, including a Madelung 
term, fixed in accordance to (Makov and Payne 1995), plus constant and a^^ fitted terms. 
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Figure 2: The multiplet spectra for ji — 2 and 3 holes, including both the e-e and the e-ph 
couplings in the anti-adiabatic approximation, as given by the total effective parameters listed in 
the last column of Table ^. The low-spin ground state for n > 3 is probably an artifact of the 
anti-adiabatic overestimation of e-ph energetics. 
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Figure 3: The n = 2 holes JT spin gap to the lowest triplet state, as a function of the dimensionless 
linear coupling strength g"^ = —4F3/{huj) of a single Hg mode (O'Brien 1972, Manini and Tosatti 
1998) (no Coulomb interaction) with pure r = 1 coupling, similar to the 271 cm~^ mode of Ceo, 
the one with the largest coupling. The same spin gap is plotted for comparison in the anti- 
adiabatic (dot-dashed - slope 5/16) and adiabatic (dashed - slope 83/320) limits. The vertical bar 
at = 12.5 locates the total effective r = 1 coupling of Ceo positive ions, according to the DFT 
estimate of Manini et al. (2001). 
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